Multipartite continuous- variable entanglement 
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Necessary and sufficient observable conditions for the nonnegativity of all partial transpositions of 
multi-mode quantum states are derived. The result is a hierarchy of inequalities for minors in terms 
of moments of the given state. Violations of any inequality is a sufficient condition for entanglement. 
Full entanglement can be certified for a manifold of multi-mode quantum states. A Mathematica 
package is given for a systematic test of the hierarchy of conditions. 
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Quantum information is a rapidly developing field of 
modern physics. Entanglement is considered to be the 
key resource for a variety of applications, such as quan- 
tum computation, quantum state teleportation and oth- 
ers. Recently increasing interest arose in the use of con- 
tinuous variable (cv) entangled states in quantum infor- 
mation processing, for a review see 0] . So far many such 
considerations deal with Gaussian entangled states, due 
to both the simplicity of characterizing their entangle- 
ment properties and the possibilities of their experimen- 
tal realization. 

Criteria of bipartite entanglement of Gaussian states 
are well established. Even necessary and sufficient con- 
ditions for the entanglement of such states arc known 
and some generalizations have been formulated 
the more general case of non-Gaussian quantum states 
a complete characterization of entanglement is still an 
open problem. A powerfull approach of characterizing 
an important class of entangled states is based on the 
Peres-Horodecki criterion which relies on the positive 
but not completely positive map of partial transposi- 
tion 0. Whenever the nonpositivity of partial trans- 
position (NPT) is identified, the state is entangled. 

Recently the necessary and sufficient conditions for the 
NPT of a bipartite quantum state have been formulated 
by using moments of arbitrarily high orders . This ap- 
proach yields a manifold of new entanglement conditions 
and unifies a series of previously known ones, including 
such whose derivation was not directly based on NPT. 
It is also of great importance that an efficient method 
has been developed for measuring the quantum correla- 
tion functions that are needed in these criteria [5(. First 
attempts have been made to extend this approach, to 
describe bipartite entanglement of more than two modes 
and entanglement beyond NPT 0, 0] • 

The characterization of multipartite entanglement is 
clearly much more sophisticated than is the bipartite 
case. For example, multipartite entanglement is con- 
sidered to be important for quantum computation. For 
some discussion of multipartite entanglement in finite di- 
mensional Hilbert spaces, cf. e.g. Q. Furthermore, mul- 
tipartite cv entanglement has been studied by using mo- 
ments of second order 0. Presently general criteria for 



multipartite cv entanglement of arbitrary quantum states 
are unknown. 

In the present contribution we derive necessary and 
sufficient conditions for the positivity of all partial trans- 
positions of a general cv multi-mode quantum state. 
Hence we obtain a hierarchy of necessary conditions for 
separability. The violation of any such condition is suf- 
ficient for entanglement. We show that full multipar- 
tite entanglement of a variety of quantum states can be 
proven by analyzing all possible partial transpositions. A 
Mathematica package for completely testing the partial 
transpositions is given, which can use the measured or 
calculated quantum correlations of a given state. 

There are several notions of separability. Let us start 
to consider full separability ^(J ■ A general n-mode state 
g is called fully separable if it is a convex combination of 
separable states, 
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where g^ is a state of the ith mode (i — 1, . . . n). The 
numbers pk satisfy the conditions pk ^ 0, J2k^°iP k = 1' 
and the series converges in the trace-norm. 

This definition can be extended in the following natural 
manner. Let tt — {I±, . . . , I p } be a decomposition of the 
set Af n into a disjoint union of the subsets Ii C fif n , i = 
1, . . . ,p: 



M n = h U . . . U I, 



v 



h n Ij = 



if ijtj. (2) 



The number of parts of the decomposition n we denote 
as \ir\ = p. The state g is called 7r-separable if it can be 
represented as 
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where g^ is a state of the multi-mode part of the orig- 
inal system formed by the modes with the indices in Ii, 
i = 1 , . . . , p, and the numbers pk satisfy the same condi- 
tions as in the case of full separability. For a related clas- 
sification of multipartite qubit states we refer the reader 
to [H|. 
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A decomposition a = { Ji, . . . , J q } is said to be finer 
than 7r, a -< n, if for any j = 1, . . . , q there is an 
i = 1, . . . ,p such that Jj C I t . The finest decomposition 
is t = {{1}, . . . , {n}}, which corresponds to the notion 
of full separability. It is clear that a 7r-separable quan- 
tum state is also tr-separable for all a such that it -< a. 
In particular, a fully separable state is 7r-separable for 
any 7r. Thus it makes sense to consider only minimal el- 
ements of the set T>(g) of decompositions n such that g 
is 7r-separable. In [12j it was shown that there are states 
g such that T)(g) has different minimal elements or, in 
other words, that some states can be separated in differ- 
ent incomparable ways. The hierarchy of the notions of 
separability considered here can be ordered as: 

full separability —* a — separability 
— > 7r — separability with er -< 7T. 

The full separability is the most restrictive one (the 
smallest number of states are fully separable). Con- 
versely, for entanglement we have: the class of entangled 
states that violate full separability is the largest, but for 
many applications not the most useful one. 

To determine whether a given multi-mode state is en- 
tangled or not is in general a highly nontrivial problem. 
A partial and powerful approach, which can verify entan- 
glement but cannot verify separability, is based on the 
notion of partial transposition. The transposition T of 
a single-mode system is defined via T(|7i)(to|) = \m)(n\, 
where {|n}} is the Fock basis in the Hilbert space asso- 
ciated with the system. For a multi-mode system one 
can introduce the notion of partial transposition, when 
only some modes are transposed and the others are left 
unchanged. 

For an n-mode system one can construct 2™ partial 
transpositions, which are in one-to-one correspondence 
with subsets / of the set Af n = {1, . . . , n}. The partial 
transposition corresponding to the subset / C J\f n we 
denote as PT/. For 1 — and / = Af n we have PT = 1 
and PT/v re = T to t respectively, where T to t is the total 
transposition of the multi-partite system. Since T 2 = 
1 we have the relation PT/ o PT,/ = PT/aj, valid for 
all /, J C J\f n , with IAJ = (I \ J) U (J \ I) being the 
symmetric difference of I and J. Taking J = Af n we get 

PT 7 = PT /0 T tot , (5) 

where / = Af n \ I is the complement of /. 

One can easily see that the partial transposition of a 
separable quantum state is again a separable state, hence 
it is nonnegative defined. One can use this fact as a test 
for entanglement: any negativity of any partial transpo- 
sition is a clear signature of entanglement. Note that the 
trivial partial transposition PT = 1 and the total trans- 
position PTjv„ = Ttot of any state are always quantum 
states, so there are 2™ — 2 choices of the subsets / C J\f n 
which can lead to negativity. Due to the relation (JSJ) and 



the fact that the total transposition of any system is al- 
ways nonnegative it is clear that for any subset / C N n 
it makes sense to test only one of the partial transpo- 
sitions PT/ and PTj, the other gives the same result. 
Finally, there are 2™" 1 — 1 partial transpositions to test, 
which can be identified with non-empty subsets of M n -i 
(including the set M n -i itself). 

It is not trivial that the study of partial transposi- 
tions gives insight in multipartite entanglement proper- 
ties. Partial transposition tests only for bipartite entan- 
glement since there are only two groups of modes: those 
we transpose and those we do not. Nevertheless, our 
approach allows one to certify full multipartite entan- 
glement of a given n-mode quantum state. Contrary to 
the notion of full separability, this requires that the state 
under study does not show 7r-separability for any decom- 
position 7r. This may be of some practical interest. In 
applications of multimode entangled quantum states the 
user may need to share n modes among k parties (fc ^ n) 
with the requirement of full entanglement. 

To certify full entanglement of some n-mode cv quan- 
tum states, it is sufficient to demonstrate that all 2™ _1 — 1 
nontrivial partial transposition are nonpositive. In fact, 
if the state is not fully entangled then there exists a de- 
composition 7r such that g is 7r-separable. According to 
relation Q it is clear that any 7r-separable state must 
fulfill at least one bipartite separability condition. Hence 
the proof that all partial transpositions are nonpositive 
is a certificate for full entanglement. The details of the 
practical realization of such a test will be given below. 

A Hermitian operator A is nonnegative if and only if 
(i I)a = /) ^ 0> f° r au operators / whose nor- 

mally ordered form exists A}. Assuming the normally or- 
dered form of / exists, one can write / = X)k ilo c ki^ k ^ 1 , 
where k = (fci, . . . , k n ), 1 = (h,...,l n ) and a^ k (a 1 ) 
means a\ kl . . . a^ fe ™ (a 1 ^ . . . a'j 1 ). The mean value (PJ)a 
is a quadratic form with respect to the coefficients Cki, 

(Pf)A= E Ck iC Pq (a tl a k atP^) A . (6) 

k,l,p,q=0 

To write this quadratic form in the standard way, let us 
order the multi-indices (k, 1) and numerate them with 
a single number. The exact order is unimportant; for 
convenience we use the following one: for any two multi- 
indices u = (k, 1) and v = (p, q) 

"<v-l' u ! < ! v ! or , (7) 

I |u| = |v| and u < v, 

where |u| = Y^t^i + Si h an d u <' v means that the 
first nonzero difference p n — k n , q n — l n , . . . , pi — ki , qi — 
l\ is positive. Combining the pairs of n-dimensional 
multi-indices (r, s) into the single 2n-dimensional ones 
as (r, s) — > (si ), we see that the order 

(|7|l becomes the graded antilexicographical order, which 
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is defined for any dimension d as follows: for two d- 
dimensional indices u = (ui, . . . , Ud) and v = (v\, . . . , Vd) 



u <s ralcx v 



|u| < |v| or 

lul = Ivl and u < alcx 



(8) 



where u < alcx v means antilexicographical order, i.e. the 
first nonzero difference Vd—Ud, ■ ■ . , v\—u\ is positive. The 
resulting ordered sequence of the moments starts as fol- 

The quadratic form 



lows: 1, (a\), (&{}, . . . , (a n ), (&1), _ 
© can now be written as (P/)a = J2t^Lo^' 



where M r 



with 



)Ai c r — Ckl, C s — upq 

(k, 1) and (p, q) being the rth and the sth indices, re- 
spectively, in the ordered sequence of moments. 

The nonnegativity of this quadratic form is equivalent 
to the nonnegativity of all its principal minors. That is, 
(/ I)a 0, if and only if the conditions 
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are satisfied for all N — 1,2, ...,oo and for all R = 
{r 1 ,...,r N } with 1 < r\ < . ■ ■ < r N . For R = 
{1, 2, . . . , In + 1} the corresponding minor is denoted by 
D2n+x and it reads as 



1 (ax) (a\) ... (at) 
(a\) (a\ax) (a\ 2 ) ... (a\at) 

(a n ) (ax(i n ) (a[a n ) . . . (a+ a„) + 1 



(10) 
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FIG. 1: Fhe minors d 1 and d?, are given for different noise 
levels: n = (solid line), n = 0.01 (dashed line) and n = 0.05 
(dot-dashed line). 



The other way around, the conditions 111|) can be formu- 
lated as follows: if there are N indices 1 ^ r\ < . . . < rjy 
and a non-empty subset / C Mn-x such that < 0, 
than the state under consideration is entangled. Note 
that, for the case of n = 2 and / = {2} the condition 
D\ < is exactly the Simon condition |2j,|^|. 
As an example let us consider the mixed state 



g(a)=N(a) / a)\^((3))W)\ d 2 /3, 



(14) 

where the pure state \ip(/3)) is defined by 

n 

|^(/3))=^|/3 1 ,...,A-i,-A,/3i+i,--^n), (15) 



For any nonnegative operator A the determinant Dm+\ 
is nonnegative: the condition D2 n +x ^ is necessary for 
nonnegativity of A (we have omitted the subscript A) . 

Now we can apply the conditions 10 to the partially 
transposed density operator A = PT/(p): it is nonnega- 
tive defined if and only if the conditions 



M 1 ... M 1 



M 1 M 1 
lyl r N r 1 r N r N 



^0 



(11) 



are satisfied for all N, where M/ t are the moments M st 
calculated on the partially transposed density operator 
PTi(g). The moments M st = (a ,1 a k a tp a q ) of any quan- 
tum state can be written as 



Mst (H";' a 



"i X 4 



(12) 



The moments M\ t of the partially transposed density op- 
erator read as 
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We/ 



a- af a- a/ 



n 



»t'« "fc "tPi "H 

a] a; a/ a? 



(13) 



and Ps(/3, a) describes noise of mean photon number n, 



Pn(M =T\Pn i (0i,a i ), Pn(P,a) = -le"^^ . 
i=i 

(16) 

This state is a noisy cv analog of the VF-state. Here 
we consider the four-mode case, n = 4. For simplic- 
ity, we may start to consider the 2x2 principal minors 

^{13,35} = ^(1,2:3,4), ^{20, 33} = d(l,3;2,4) an d d{22, 31} = 

d(2,3 ; i,4), where duj-k i) is defined via 



(al&iajaj) (d\d]d k di) 
(aidja\a\) (d\dkd\di) 



(17) 



This principal minor corresponds to the operator / = 
cxdidj +C2&kai in the expression JBJl. The behavior of the 
partially transposed minors dx — df^ 3g , = d\ 2 ^ 

{3} _ ,{1,2,3} , , _ ,{1,2} _ ,{1,3} 



'{13,35} ~~ 

JT.*r _ j{1,2,3} , , _ ,{1,2} _ ,{1,3} _ ,{2,3} 
a { 13,35} — a {13,35} ^ nCL 0,2 ~ fl {13,35} — "{20,33} — ^{22, 31} 

Note that these minors coincide 



dYSL, = dYSL, = d) 

is shown in Fig. 
only for our highly symmetric case (a, = a and nt = n 
for i = 1,2,3,4). Figure n illustrates the behavior of 
the minors dx and c?2 as function of \a\ for different noise 
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levels. It shows that the maximal negativity can be found 
for the noiseless case (n = 0). One can see that for 
a range of \a\ (which depends on the noise) the state 
under study is fully entangled, what can be demonstrated 
already by lower-order moments. For larger values of \a\ 
or for a larger noise level lower-order moments do not 
show the negativity. 

To certify full multipartite entanglement, one needs a 
systematic procedure to deal with the data of a given 
quantum state to be analyzed. The data may be cal- 
culated or measured moments, depending on whether a 
theoretical or an experimentally measured quantum state 
is considered. Let us present in detail the algorithm for 
computing the nth multi-index in the graded antilexico- 
graphical order. First, for any d-dimensional index u let 
us find the next multi-index v, which has the number 
n + 1, i.e. such an index v that u <s ralex v but there 
are no other indices w with u <s ralox w <;graicx v rp nc 
index v, next to the given u, can be computed with the 
following algorithm: 
nextMultiindex(u) 
input: u = (m, . . . ,u d ) 
output: minimal v such that u <s ralcx v 

1. find the minimal i such that ui ^ 0; if all it, = 
set i = d. 

2. if i = d 

v= (u d + l,0,...,0); 
else 

V = (Ui - 1, 0, . . . , 0, u l+1 + 1, u l+2 , . . . , u d ). 

" * ' 

i 

The nth multi-index can be computed using this func- 
tion repeatedly starting with the minimal multi-index 
(0,...,0): 

nthMultiindex(d, n) = nextMultiindex^ 11 (0, . . . , 0), 

d 

(18) 

where / [OI (u) = u and /M(u) = /(/ [ ™" 1] (u)) for all 
n > 0. In the case of d = 2 one can get the explicit 
expression for the nth two-dimensional index: 

nthMultiindex(2, n) = 

f(N + l)(N + 2) N(N + 1) \ (19) 

n, n 1 , 

\ 2 2 J' 

where N = \ y/Sn + 1 ~ 3 ], with \x~\ being the smallest inte- 
ger greater than or equal to x. 

Let us comment on the realization of the measurements 
to verify full multipartite entanglement. Our chosen ex- 
ample of a fully entangled quantum state clearly shows 
the strong sensitivity of entanglement with respect to 
noise effects. Any quantum-state reconstruction based 
on homodyne detection with imperfect detectors includes 
Gaussian noise effects, for a review see 13] , that may pre- 
vent one from verifying full entanglement. Therefore we 



have recently proposed the method of balanced homo- 
dyne correlation measurements to overcome such prob- 
lems Upon properly balancing the setup, any kind 
of minor to be considered in our method is simply pro- 
portional to a corresponding power of the detection ef- 
ficiency. Though different moments are proportional to 
different powers of the quantum efficiency, in the minors 
they are combined in such a way that the whole minor 
is proportional to one common power of the efficiency. 
Hence all the inequalities to be used to verify the non- 
positivity of partial transpositions are only multiplied by 
a positive factor, so that the method also works in the 
case of imperfect detection. 

In conclusion, we have applied the partial transposition 
approach to get insight into multipartite entanglement of 
mixed multi-mode cv quantum states. The method al- 
lows us to certify full entanglement of a broad class of 
multimode quantum states, also of non-Gaussian ones, 
which requires making use of higher-order moments. We 
also provide a Mathematica package to automate the pro- 
cess of testing partial transpositions, given that the nor- 
mally ordered moments are defined either analytically or 
as a result of measurements. 
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